Nontrivial temporal scaling in a Galilean stick-slip dynamics 
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We examine the stick-slip fluctuating response of a rough massive nonrotating cylinder moving 
on a rough inclined groove which is submitted to weak external perturbations and which is main- 
tained well below the angle of repose. The experiments presented here, which are reminiscent of 
the Galileo's works with rolling objects on inclines, have brought in the last years important new 
insights into the friction between surfaces in relative motion and are of relevance for earthquakes, 
differing from classical block-spring models by the mechanism of energy input in the system. Robust 
nontrivial temporal scaling laws appearing in the dynamics of this system are reported, and it is 
shown that the time-support where dissipation occurs approaches a statistical fractal set with a 
fixed value of dimension. The distribution of periods of inactivity in the intermittent motion of the 
cylinder is also studied and found to be closely related to the lacunarity of a random version of the 
classic triadic Cantor set on the line. 

PACS numbers: 81.40.Pq, 05.40.-a, 05.45.Df, 05.70.Np 



I 

C 

O 
o 



(N 
> 
m 
oo 

O 
in 
o 

c3 



i 

-a 
c 

o 
o 



X 



I. INTRODUCTION 



It can be conjectured that the experimental verifica- 
tion of fluctuations in the dynamic behavior of a block 
on an incline was first verified in the beginning of the 
17th century. It is known that as early as 1604 Galileo 
believed in the now classical scaling law s ~ t 2 connecting 
the distance s and the time t for an object falling close to 
Earth under the action of gravity Q. Galileo's applica- 
tion of this law for rolling objects on inclined planes was 
discussed in his last book The Two New Sciences, pub- 
lished in 1638. However, no explicit or implicit mention 
of the application of this simple scaling law for nonrotat- 
ing objects is made in Galileo's works. Certainly Galileo 
tried to study the pure translational motion of objects 
on an incline, but he was possibly unable to encapsu- 
late the experimental results in a simple mathematical 
law as the previous one: the experimental outcomes in 
this case are very complex, apparently nonreproducible, 
due to the fluctuations in the friction force acting at the 
incline— sliding object interface. Few nonspecialized text- 
books have called attention to the inherent fluctuating 
character of friction forces between solid surfaces (an ex- 
ception is Ref. ;2|). Only recently friction fluctuations 
involving solid surfaces have been quantitatively studied 
in a number of macroscopic situations 0- 01 0] ■ 

In the last few years, it has been shown that inter- 
mittent sliding or stick-slip dynamics of a rough solid 
nonrotating cylinder on a rough inclined groove submit- 
ted to small controlled perturbations is a fluctuation phe- 
nomenon characterized by nontrivial spatiotemporal scal- 
ing laws |6J, [jj and complex critical exponents |8| if the 
inclination is well below the angle of repose. In partic- 
ular, the time series of intermittent slidings associated 
with the stick-slip motion of the cylinder on the incline 
were found to present many similarities to time series of 
earthquakes: The sliding distribution is described by the 



Gutenberg and Richter law n(s) ~ s , where n(s) is 
the number of events of size s, and the Omori law for 
the number of smaller events occurring at a time t after 
a large event, n(t) ~ t~ p , where the exponent p is an 
anomalous one, lying between 0.25 and 0.45, and may 
be a complex number. Stick-slip dynamics is a topic of 
broad interest, and it traditionally refers to the situa- 
tion in which a solid on a horizontal surface is pulled 
at a constant driving velocity 0. In the experiments 
discussed in this paper, the stick-slip dynamics appears 
as a consequence of a completely different mechanism: to 
start the slip we resort to small mechanical perturbations 
on an inclined surface, whose angle with the horizontal 
is below the angle of repose. It is important to notice 
that this mechanism constitutes indeed the basic differ- 
ence from many models of earthquakes, which use blocks 
and springs to simulate the motion of tectonic plates 
[ToL [Til ITsj . A solid body on a perturbed incline is an 
example of a nonequilibrium system receiving an incom- 
ing energy flux. If energy is continuously injected into 
nonequilibrium systems, a complex sequential response 
characterized by time series of events of all sizes is of- 
ten observed. Besides sliding blocks on inclines, other 
examples of systems and phenomena associated with a 
similar type of temporal fluctuati ng r esponse are piles of 
sand and other granular materials |l3lll4| : acoustic emis- 
sion from volcanic rocks and microfracturing processes in 
general [HI Il6[ ; interface depinning in magnetic systems 
|17|: stick-slip motion in lubricated systems and tur- 
bulence 0| , among others. 

This work reports on basic geometric aspects of the 
fluctuating behavior observed in experimental time series 
of sliding events of a rough nonrotating metallic cylin- 
der moving intermittently on a rough groove weakly per- 
turbed by external impacts and maintained well below 
the angle of repose. The Gutenberg & Richter law dis- 
cussed for this system in Ref. is a relation between 
two macroscopic variables, independent of the sequence 
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of the sliding events; furthermore, the Omori law studied 
in Ref. Q and the Hurst exponent analysis presented in 
Ref. are defined for a mesoscopic scale associated with 
time windows of many slips within a sequence of sliding 
events. In this work, we investigate statistical variables 
which are strongly related to the microscopic details of 
the time series of slidings. The experiments giving ori- 
gin to these time series are very time consuming but the 
results are quite robust in the sense that they are inde- 
pendent of the particular value of several experimental 
parameters. In particular, we learn from our experiment 
that (i) robust and nontrivial temporal scaling laws ap- 
pear in this stick-slip dynamics, and (ii) the time-support 
where dissipation occurs is also statistically robust, defin- 
ing a set whose dimension lies in the interval 0.61 ± 0.01. 
Moreover, (iii) the lacunarity function 19] describing the 
statistical distribution of time domains free of dissipation 
is shown to be related to the corresponding lacunarity of 
a random version of the classic triadic Cantor set on the 
line. The Galilean experiments giving rise to the inter- 
mittent time series studied in this work are described 
in the next section. In Sec. 3 our results are presented 
and discussed. We then justify the use of mathemati- 
cal constructions as the Cantor sets on the line to give a 
statistical description of the intermittent distribution of 
slidings on the time axis. Conclusions are made in Sec. 4. 



II. EXPERIMENTAL DETAILS 

The basic apparatus used to obtain the experimental 
data discussed in this work consists of a rigid V-shaped 
anodized aluminum chute made of a corner plate of 5 mm 
thickness and with a 90° angular aperture symmetrically 
disposed with respect to the vertical plane (Fig.^l. The 
chute is rigidly maintained with an inclination 9 with 
respect to the horizontal, and it is supplied with an ar- 
ticulated hammer of mass m which hits the base of the 
chute with a controlled (fixed) velocity. The results re- 
ported here are statistically invariant for m in the interval 
studied from 75 g to 175 g. The experimental apparatus 
relaxes elastically within the time of 5 s between two con- 
secutive impacts of the hammer. The system is mounted 
on a 200 kg table isolated from mechanical vibrations. 
On the chute (with an effective length Lq = 1500 mm) is 
placed a metallic cylinder of length L (5—1000 mm), and 
the system chute + cylinder operates significantly below 
the critical angle of repose 9 C — tan _1 /x s , where fj, s is the 
coefficient of static friction of the cylinder on the chute. 
The weak impacts of the hammer do not lead to jumps of 
the cylinder out of contact. In all experiments, the incli- 
nation was in the interval 12°-18°, with 16° < 9 C < 32°, 
and the reduced angle (9 C — 9)/9 c was typically in the 
interval 0.28-0.35. 

The initial condition in all experiments refers to the 
cylinder at the top of the groove. Afterwards, the num- 
ber of the impact of the hammer after which each sliding 
occurred, i, is recorded, as well as the corresponding slid- 
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FIG. 1: Schematic diagram of the experimental apparatus: 
the incline A has a length Lq = 1500 mm, and is rigidly 
maintained at an inclination 8 < 9 C (=angle of repose) with 
the horizontal. The cross section of A is shown in B; C is 
a cylinder with length L varying from L — 5 to L = 1000 
mm. Sliding events occur after a variable number of impacts 
of the hammer D. The time series discussed in the text and 
illustrated in Fig. [5] are chains of 100 sliding events. 



ing length A(t). In the regime of inclination used in the 
experiments, the series of induced sliding events are in- 
termittent, i.e., a fluctuating number of many controlled 
perturbations of the hammer (time units) is necessary to 
induce a single sliding event of the cylinder. The total 
duration T of each time series corresponds to the number 
of hammer impacts after 100 sliding events A ^ 0, and 
varies from 364 to 6, 979 time units or perturbations [20| . 
The number of sliding events was fixed in such a manner 
that in all experiments the time series were recorded be- 
fore the cylinder hit the end of the incline. Thus, each 
sliding sequence refers to a single sweep of the cylinder 
along the groove. The total number of perturbations in 
our experiments exceeded 56, 000, corresponding to 6, 000 
sliding events. For illustration, Fig. [2] shows typical in- 
termittent sliding series for two values of the length L of 
a massive aluminum cylinder of diameter 9.6 mm. 

III. RESULTS AND DISCUSSION 

The dynamics of translation of a solid body on an in- 
clined groove is dependent on the angle of inclination 9: 
(i) If 9 > 9 C , there is a trivial single continuous sliding 
event (with or without perturbation on the groove), (ii) 
If 9 is slightly smaller than 9 C , the continuous phase (i) 
is replaced by a discrete homogeneous phase with a one- 
to-one correspondence between perturbation and sliding 
response 0. In this case, there are spatial fluctuations in 
the magnitude of the sliding, but no temporal gaps be- 
tween sliding events. This phase can be named Euclidean 
because the support of dissipation has the same dimen- 
sion D = 1 of the time axis, (iii) If 9 is significantly 
smaller than 9 C (as quantified in the previous section), 
there are both spatial fluctuations in the magnitude of 
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FIG. 2: Typical nonhomogeneous intermittent time series of 
sliding events for massive cylinders of aluminum of length (a) 
L = 20 mm and (b) L = 200 mm, with (0 C - 6)/6 c w 0.30. 



the sliding events and fluctuations in the time between 
consecutive sliding events, (iv) For very small angles, 
the resting state of the cylinder is obtained. It is the 
nonhomogeneous intermittent regime (iii) that is studied 
in the present work: Fig. [5] shows the nonhomogeneous 
intermittence observed in two time series obtained from 
our experiments. This is interesting and nontrivial be- 
cause in principle the dissipation could present gaps ho- 
mogeneously distributed along the time. However, this 
intermittent homogeneous-gap phase was never detected, 
i.e., nature seems to prefer the heterogeneous clustering 
of sliding events. Perhaps, this homogeneous gap phase 
could be detected in experiments using a single crystal 
placed on another single crystal. 

In an experiment like the block— incline studied here, 
where many variables and parameters are involved, it is 
of crucial relevance to know which elements related to 
the apparatus are of importance for the behavior of the 
system, and how they could influence any finding. For 
example, we know that cylinders of smaller length take 
more time (more impacts of the hammer are needed) to 
slip, and this reflects the fractal nature of the rough- 
ness of the surfaces in relative motion ppj | . Furthermore, 
from a "first sight" into our experiment, one could be 
easily led to the conclusion that the sliding susceptibil- 
ity of the blocks should increase in time since they get 
closer to the impact source (the hammer) after each slip 
(see Fig. [Q. However, figure [3] shows that this is not 
the case. In the main plot of this figure, we show the 

cumulative sum A t — X/«=i •M'O °f the magnitude A of 
the sliding events as a function of time t for a cylinder of 
length L = 10 mm, while in the inset the same quantity 
is shown for a cylinder of length L = 100 mm. As we 
can see from these plots, A t increases quite linearly in 
i, i.e., the sliding probability remains constant in time, 
and this result is independent of the length of the cylin- 
der. In other words, the energy released by the hammer 
impact at the bottom of the apparatus is transmitted in 
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FIG. 3: Cumulative sum A t = Yli=i M^) °f the sliding events 
as a function of time for two different experiments: In the 
main plot, the corresponding time series was obtained with a 
cylinder of length L = 10 mm, while in the inset the length 
of the cylinder used in the experiment was L = 100 mm. The 
dashed lines in the two plots are guide to the eye meant to 
show the linear behavior of At. 



average homogeneously along the chute. Figure |31 has to 
be distinguished from Figs. 1 — 4 of Ref. where we 
plot the number of sliding events as a function of time, 
despite their magnitude. On the other side, the behavior 
of the system may change dramatically depending for in- 
stance on the values of the inclination of the chute and on 
the intensity of the hammer impact. Decreasing [increas- 
ing] the angle 8 of the incline and/or the intensity of the 
perturbations to values far from the ranges mentioned in 
the last section would bring the system from the regime 
(iii) defined in the last paragraph to the regime (iv) [to 
regime (ii) or even (i)]. However, in this paper, we are 
just interested in the regime (iii) of nonhomogeneous in- 
termittent motion of the cylinder on the chute. In the 
next paragraphs we will quantify this intermittent dis- 
tribution by calculating the dimension of the dissipation 
support and its lacunarity. 

After assigning "1" ["0"] to each time position in which 
the cylinder slides [is at rest] , we define the temporal dis- 
sipation support £ as the subset of the time axis consist- 
ing of l's. To quantify the time distribution of the sliding 
activity, we count the number N(t) of one-dimensional 
boxes of size r necessary to cover S. Here, the vari- 
able t is defined as At/T, where At is the corresponding 
length of the box in time units, and T is the total num- 
ber of perturbations. In the inset of Fig. 0] we show 
N(t) versus r for some typical series. The circles shown 
in the main plot represent the complete ensemble aver- 
age (N(t)), which was obtained by dividing the r axes 
into 15 equal intervals and calculating the mean values 
of N(t). 
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FIG. 4: Log-log plot of the average number (N(t)) of one- 
dimensional boxes of size r necessary to cover the regions 
along the time axis where the sliding activity is concentrated. 
The average (•} is over the entire experimental set of time se- 
ries of sliding events (circles). The asterisks denote the same 
quantity for an ensemble of equal size of synthetic series: ran- 
dom Cantor sets whose number of events is equal to the num- 
ber (100) of sliding events (see text). The dotted line has 
as slope the golden mean (V5- l)/2 ~ 0.618.... The inset 
shows N(t) for typical experimental time series. 



For comparison, we show with asterisks in Fig. 0] the 
corresponding function (Nurcs(t)} averaged over an en- 
semble of different random Cantor sets defined as fol- 
lows. First, the algorithm for the triadic random Can- 
tor set described in Ref. [2lJ is implemented for a line 
with length T = 3 7 = 2187 (this value is close to the 
mean size of the dissipation support of the series stud- 
ied, (T) = 2175): each line segment is divided into three 
equal parts and one of them is thereafter randomly re- 
moved. After the 7th iteration, this random construc- 
tion contains just 2 7 = 128 events. We then take out 
28 events randomly, to get a modified random Cantor 
set (MRCS) with the same number (100) of events as in 
the experimental series. (Nmrcs(t)) scales as t ~ Dmrcs , 
with Dmrcs = 0.61 ± 0.01, as shown by the dotted line 
in Fig. The continuous line in the plot has the same 
slope as the dotted one. As we can see, this line gives 
a good approximation to the corresponding slope of the 
experimental trend line along almost two decades in r. 
Although the experimental data clearly do not present 
a perfect scaling - due to the low dimensionality of the 
phenomenon and the consequent difficulty to obtain se- 
ries with many events — (N(t)} indicates within the error 
bars that the dissipation support approaches a set with 
dimension Dmrcs = 0.61 ± 0.01, which is the fractal di- 
mension of the random Cantor set defined above. It is im- 
portant to notice that the scaling relation shown in Fig.0] 
does not mean that the dissipation support is a fractal 
set. Moreover, all the random Cantor sets studied here 



are statistical fractal sets and have in average the same 
finite size (T) of the experimental series, as well as the 
same number of events as in the corresponding sliding se- 
quences. A quite general random Cantor set construction 
was introduced by Mauldin and Williams |22j . It consists 
of the recursive division of the line segment into three 
parts of different sizes, which are determined by two ran- 
dom parameters < r%,r2 < 1, followed by the removal 
of the central segment. The dimension of the random 
Cantor sets obtained with this algorithm is known as the 
golden mean, Dmw = {V§ — l)/2 = 0.618 . . .. Mauldin 
and Williams showed that although the infinite combina- 
tions of ri and r2 lead to completely different geometries, 
the resulting random Cantor sets appear robust with re- 
spect to their dimension. The value of -Dmw is universal 
and close to the value -Dmrcs = 0.61±0.01 obtained from 
Fig. 0] We speculate that this general class of random 
recursive dynamics introduced by these authors could be 
associated with the physical phenomenon discussed in the 
present work. 

An additional method to define a fractal dimension is 
based on the mass-size relation 23]. In our case this 
means to count the number of sliding events, M(r), 
within an interval of time r. For a scaling distribution of 
intermittent events on the time axis we expect that the 
ensemble average (M(r)) ~ r* 5 , where the mass-exponent 
5 < 1. In Fig. Owe show that, in fact, (M(r)) behaves as 
a power law of r with the exponent <5 = 0.65±0.05, which 
is essentially the same as -Dmrcs within the fluctuation 
bars. Again, this figure has to be distinguished from 
Fig. 3 in Ref. where we show the number of events 
between two large slips, which provides a mesoscopic de- 
scription of the time series, as opposed to the variable 
(M(t)), which quantifies the microscopic distribution of 
sliding events along the dissipation support. 

Another way to get randomized Cantor sets in our 
problem would be taking out randomly 28 of the 128 
events of a seventh order triadic Cantor set (TCS), whose 
fractal dimension is -Dtcs = In2/ln3 = 0.6309... (ob- 
tained by recursive removal of the central third of each 
segment 0,|23|). Let us define a modified triadic Can- 
tor set (MTCS) of seventh order as a triadic Cantor set 
after the random removal of 28 events. It can be easily 
verified that this operation leads to a set with dimension 
£> MTCS = 0.61 ± 0.01. This means that both the MRCS 
and the MTCS simulate, from the point of view of fractal 
dimension, the dissipation support of the time series of 
sliding events. However, it is clear that these two kinds 
of Cantor sets are very different in respect to their ori- 
gin. Clearly, the MRCS is more random than the MTCS. 
To remove the degeneracy between these two sets due to 
their common value of fractal dimension, we examine an- 
other statistical property: the lacunarity function. 

Lacunarity, a concept introduced by Mandelbrot, tries 
to quantify the texture of a fractal (23 ■ Here we use the 
mathematical algorithm introduced by Allain and Cloitrc 
[T^ | to calculate the lacunarity of the dissipation support 
of the sliding events. A short description of their method 
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FIG. 5: Log-log plot of the average number (M(t)) of sliding 
events within an interval of time r. The average (.} is over the 
entire experimental set of time series of sliding events (circles). 
The line has slope S = 0.65 ± 0.05, which is interpreted here 
as the (mass) dimension of the statistic fractal set defined by 
the dissipation support of the experimental time series. This 
value is in agreement, within the fluctuation bars, with the 
dimension found in Fig. 0] calculated with the box counting 
for the same values of the time interval r. 
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FIG. 6: Log-log plot of the lacunarity function (A(r)) for the 
experimental time series (circles). The dashed line gives the 
same function for an ensemble of MTCS of order 7 (2 7 events, 
T — 3 7 ) with 28 events randomly eliminated (see text). The 
inset shows A(r) for some typical time series. (A(r)} decays 
effectively as the power law r ~ ' 36 for 3 x 10~ 3 < r < 2 x lO^ 1 
(continuous line). 



can be given as follows. First, a time window of size r 
is denned, which is translated by R(r) = T(l — r) + 1 
steps (time units), from the beginning up to the end of 
the time series, which has size T. The number n(s, t) is 
defined as the number of times in which the window of 
size r contains s sliding events during translation. Thus 
the probability that the window encounters s occupied 
time units during translation is p(s,r) — n(s,r) / R(r). 
It follows that the lacunarity A(r) can be defined as 

A(r) = ^ s2 P^ T \ . (l) 

EI=i s p( s > r )] 

This function has very interesting properties, being equal 
to unity when r = 1 . The inset of Fig. [5] shows the la- 
cunarity A(t) calculated for some time series studied. In 
the main plot we show the experimental ensemble aver- 
age value (A(t)), represented by circles, together with 
the corresponding fluctuation bars, as a function of r. 
As we can see from this figure, the dashed line follows 
very closely the trend of the curve defined by the circles. 
This line is no best fit obtained from the experimental 
data: it is the lacunarity calculated for the MTCS of 
seventh order of the previous paragraph. The agreement 
of the curves is noticeable. It can be seen from Fig. H3 
that (A(t)) decays effectively as the power law r~ 36 , 
along approximately two decades of variability in time, 
for 3 x ICT 3 < r < 2 x 1CT 1 , as indicated by the contin- 
uous line. 

The Cantor sets discussed here are mathematical struc- 
tures with topological dimension d t — 0, fractal dimen- 



sion D in the interval (0, 1), and embedding dimension 
d = 1 |2^|. As shown in Fig. |2 intermittence does 
not distribute homogeneously on the time axis, i.e., its 
temporal support may, in principle, be characterized by 
a dimension D < 1. Thus Cantor sets are natural 
candidates to model nonhomogeneous intermittence and 
other fluctuation phenomena defined along the ID time 
axis. Brownian motion is the prototype of the fluctua- 
tion phenomenon: if we plot the position X of a parti- 
cle executing Brownian motion in one-dimensional (ID) 
as a function of the time t, we obtain a fractal record 
X(t) having dimension Db = 3/2 on the position-time 
plane. The zeroset for the Brown function X(t), which 
is defined as those instants t for which X(t) cuts the 
time axis, is also a Cantor set, but one with dimension 
D = D B - 1 = (3/2) - 1 = 0.5 23]. Analogously, the 
temporal dissipation support studied in the present work 
can be seen as a zeroset associated with the sliding length 
distribution X(t), albeit a zeroset with larger (and equally 
robust, as in the Brownian case) dimension in the interval 
0.61-0.65. 

Finally, we must ask why indeed should this type of 
experiment exhibit nonhomogeneous intermittence with 
a random response of sliding events? A possible explana- 
tion is that the number and the individual characteristics 
of the several contacts at the cylinder— chute interface 
fluctuate continuously on the time scale associated with 
the hammer impacts. Eventually, e.g., a single interface 
contact providing a force equilibrating a large (small) 
fraction of the component of the cylinder's weight along 
the chute fails and, as a consequence, a corresponding 
large (small) resulting force appears, with the cylinder 
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tending to stop after a long (short) sliding length. The 
motion of the cylinder is followed by the rest state when 
the fluctuating interfacial roughness is able to provide 
a pattern of contacts producing the necessary opposing 
force to equilibrate the weight's component down the 
chute. Evidently as the angle (or the component of 
the cylinder's weight down the chute) decreases, it be- 
comes on average increasingly difficult to modify the rest 
state of the cylinder; i.e., the sliding susceptibility de- 
creases [2(j ■ The observed intermittence with fluctuating 
gaps denoting the interval between sliding events (dissi- 
pation) would be a manifestation of the fluctuation in 
the ability of the cylinder— chute contacts in maintaining 
the equilibrium of forces for the cylinder. This situation 
is reminiscent of the fluctuation-dissipation theorem: If 
the contacts and the corresponding forces at the interface 
cylinder— chute do not fluctuate, the cylinder remains 
at rest, i.e., there is no dissipation or sliding event; on 
the contrary, for increasing fluctuations at the interfacial 
forces, the dissipation also increases, i.e., a cascade of 
intermittent sliding events is obtained as the associated 
response. 

IV. CONCLUSION 

We found robust nontrivial temporal scaling laws re- 
lated with the complex time-distribution of sliding events 
observed in an extensive experimental study of Galilean 
stick-slip dynamics on rough inclined surfaces. It is found 
that the time-support where dissipation occurs in these 
processes approaches a statistical fractal set character- 
ized by a fixed dimension, D = 0.61±0.01. The lacunar- 
ity function for the time series of events is closely related 
to the corresponding lacunarity of a modified random 
version of the classic triadic Cantor set [2j| . It has to be 
emphasized that the geometric and statistic properties 
exhibited by the time series of sliding events are a conse- 
quence of the particular nature of the microscopic friction 
at the block— incline interface. No spatial and tempo- 
ral correlations as shown here and as those reported in 
|E SB 13 are found in these series if the block is moved 
back to the top of the incline after each sliding event, 
which would mean removing the strong correlation be- 



tween the configuration of the contacts at the interface 
and the past history of the block on the incline. Due 
to the very time-consuming nature of the experiments 
described in this paper, many aspects remain to be in- 
vestigated in the future. It is important to notice that 
the results presented in this paper are restricted to time 
series obtained from the intermittent stick-slip motion 
of the cylinder, which means the regime (iii) defined in 
section 3, and thus constrains the values of experimen- 
tal parameters as, for instance, the inclination 6 of the 
chute and the mass of the hammer to the ones mentioned 
in section 2, as said before. However, new experiments 
exploring systematically the effect of the intensity of the 
perturbation on the groove, as well as the variation of the 
materials involved and the effect of the geometry of the 
sliding objects are necessary for a better understanding 
of this fluctuation phenomenon. Furthermore, another 
important aspect to be considered is the possibility of a 
universal value D = 0.61 ± 0.01 for the dimension of the 
dissipation support in any experiment involving sliding of 
rough solids of the type discussed in this paper. In a re- 
cent paper j24j , Feng and Seto have analyzed several time 
series of acoustic emissions obtained from microfractur- 
ing in rocks, which is a process intimately connected with 
the dynamics of earthquakes, and have found that trig- 
gering events associated with a failure define a dissipation 
support with values of dimension close to the exponents 
found for our time series of slidings. The application of 
the analysis reported in the present work to the temporal 
response of other nonequilibrium systems as those men- 
tioned in the end of the second paragraph could clarify 
important aspects of universality in nonequilibrium dis- 
sipative systems. 
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